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ABSTRACT 

Powerful outflows from star-forming galaxies are expected to push away the neutral intergalactic 
medium (IGM) around those galaxies, and produce absorption- free gaps in the Lya forest. We analyze 
the abundance of gaps of various sizes in three high resolution spectra of quasars at z ~ 3 — 3.5. 
The gap statistics agrees well with a model in which galactic halos above a minimum mass scale of 
-W m i n ~ 10 10 Mq produce bubbles with a characteristic radius of Rb ~ 0.48/i _1 Mpc. Both numbers are 
consistent with naive theoretical expectations, where the minimum galaxy mass reflects the threshold for 
infall of gas out of a photo-ionized IGM. The observed gaps are typically bounded by deep absorption 
features as expected from the accumulation of swept- up gas on the bubble walls. 

Subject headings: large-scale structure of the universe — intergalactic medium — quasars: absorption 
lines — cosmology: theory 
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1. INTRODUCTION 

The formation and evolution of galaxies is regulated by 
their feedback on the surrounding intergalactic medium 
(IGM) . Radiative feedback owing to photo-ionization heat- 
ing by a cosmic UV background or hydrodynamic feedback 
owing to outflows driven by supernovae or quasars, can 
suppress the infall of IGM gas onto low-mass dark-matter 
halos (see, e.g., Ikeuchi 1986; Rees 1986; Babul & Rees 
1992; Efstathiou 1992; Thoul & Weinberg 1996; Scanna- 
pieco & Oh 2004). The feedback introduces a minimum 
halo mass above which halos can efficiently accrete gas 
from the IGM and host star formation or quasar activity. 
Infall of gas is suppressed if the gravitational potential well 
of a halo is shallow, explaining why dwarf galaxies have a 
much lower abundance than expected for low-mass halos 
(e.g., White & Frenk 1991; Kauffmann, White, & Guider- 
doni 1993). 

In this Letter, we examine a simple observational 
method to reveal the scars left by hydrodynamic out- 
flows from galaxies on the IGM. Powerful outflows from 
galaxies are expected to push away the neutral intergalac- 
tic medium (IGM) around those galaxies, and produce 
absorption-free gaps in the Lya forest towards a back- 
ground quasar (Theuns et al. 2001). Studies of the Lya 
forest — the narrow absorption lines produced by inter- 
galactic neutral hydrogen in quasar spectra — had been 
very successful in mapping large scale structure in the uni- 
verse at high redshift (for a review, see Rauch 1998). While 
the IGM dynamics is shaped by gravity alone on large 
scales (see, e.g., Cen et al. 1994; Hernquist et al. 1996; 
Meiksin & White 2001; McDonald et al. 2004), outflows 
can play an important role on small scales by introduc- 
ing gaps in the Lya forest that can be identified in high 
resolution observations of quasar spectra. 

Our model is simple: an outflow from a galaxy evacuates 
a bubble in the IGM, leaving behind little or no neutral 
hydrogen as the swept-away gas piles up on the bubble 



surface. Such a bubble would appear as a gap in the Lya 
forest absorption spectrum of a background quasar, i.e., 
a flat spectral segment with little or no absorption. The 
expected distribution of absorption-free gaps can be calcu- 
lated from the number density of halos above the minimum 
mass, as only those halos are capable of hosting supernova 
or quasar-driven outflows. For the halo mass range of in- 
terest, ~ 1O 9 -1O 11 M , the bubble radius is expected to 
be almost independent of halo mass at a given redshift 
(Furlanetto & Loeb 2003). Hence, our analytic model has 
only two free parameters: the minimum halo mass for the 
production of outflows, M m j n , and the characteristic bub- 
ble radius Rb for halos above that mass. 

There have already been some observational reports 
about the signatures of outflows around massive starburst 
galaxies (with M » M min ) at high redshifts. Shapley 
et al. (2003) discovered large scale outflows with veloci- 
ties up to ~ 600 km s -1 by studying ~ 10 3 Lyman Break 
Galaxies (LBGs) at z ~ 3. Using a large sample of quasar 
spectra which pass by foreground LBGs, Adelberger et 
al. (2003) identified the lack of neutral hydrogen within 
~ 0.5/i _1 comoving Mpc around LBGs. Recent observa- 
tions of high redshift, submillimeter galaxies also indicated 
the existence of large scale outflows (see, Genzel et al. 2004 
and reference therein). The properties of active bubbles 
around starburst galaxies may be different than those of 
mature bubbles around more common galaxies for which 
the wind material had sufficient time to accumulate on the 
bubble walls and decelerate after evacuating the bubble in- 
terior. In simulations, Croft et al. (2002) and Kollmeier 
et al. (2003) also found absorption-free gaps in their sim- 
ulated Lya spectrum caused by galactic winds. Our goal 
in this Letter is to characterize the properties of the typ- 
ical relic bubbles around the most common (M ~ M m i n ) 
galaxies at redshift z ~ 3. 

Our paper is organized as follows. In §2 we briefly de- 
scribe our analytic model and derive its basic equations. 
In §3 we compare the model with three high-resolution 
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quasar spectra. Finally, we conclude with a discussion 
about the assumptions of our model and the implications 
of our results. Throughout the paper, we assume the 
concordance cosmological model with 17m = 0.27, Cl\ = 
0.73, fi fc = 0.044, ft = 0.71, andcr 8 = 0.84 determined by 
recent WMAP data (Spergel et al. 2003). 

2. MODEL 

In our simple model, a galactic wind sweeps-up the sur- 
rounding intergalactic medium and eliminates HI absorp- 
tion within a bubble radius Rb of order a few hundred kpc 
(Adelberger et al. 2003; Furlanctto & Locb 2003). The 
bubble spends most of the time approaching a fixed co- 
moving distance from its host galaxy, at which the (decel- 
erating) expansion speed asymptotes to the local Hubble 
velocity (see section SS4). Such a bubble would produce 
a gap in the Lya forest imprinted on the spectrum of a 
background quasar. The absorption-free gap should have a 
high value of the mean transmission F [defined as exp (— t) 
where r is the optical depth] that is close to unity. Our 
goal is to calculate the expected size distribution of such 
gaps in the Lya forest and compare it with observations. 

Given a comoving density of halos per unit mass of 
dnh{M)/dM and the assumption that the wind from each 
halo of mass M produces bubble of a comoving radius of 
Rb(M), the probability of finding an absorption- free gap 
(with F ~ 1) of comoving size in the range between L and 
L + dL per infinitesimal comoving path- length, dl, is given 

by 

dP f drih dcr , . 

dLdl = J Mm . n dMdL ' ( ' 

Here, a — irb 2 is the cross-section of the bubble for an im- 
pact parameter b — \J R% — (L/2) 2 , and M min is the min- 
imum halo mass into which IGM gas accretes efficiently. 
The differential cross-section for having a bubble segment 
of length between L and L + dL is then, 

da _ f irL/2 0<L< 2R b (M); 

dL ~ \ L > 2R b (M) W 

The expected number of absorption-free gaps with a size 
larger than L is 

dP f°° r 2R b(M) j j 

dl JM min Jl dM dL' 

-£/(«-t)&- <»> 

Since Rb{M) has a rather weak dependence on M across 
the relevant range of halo masses (see Fig. 1 in Furlan- 
etto & Loeb (2003)), we assume for simplicity a single 
value for Rb when fitting observational data at a given 
rcdshift. We adopt the Press & Schechtcr (1974) mass 
function of halos and use M min as the second adjustable 
parameter. The concordance cosmological model yields 
the following values for the cumulative number density: 
n h (> M min ) = 1,0.1, and 0.01 h 3 comoving Mpc~ 3 , for 
M min = 1.3 x 10 9 , 1.2 x 10 10 , and lO n M , respectively. 
Figure 1 shows results from equation (3) for the three 
cases of M m i n = 1.3 x 10 9 M Q (dashed line), M m i n = 
1.2 x 1O 1O M (dotted line), and M min = 10 11 M (solid 
line). 



3. OBSERVATIONS 
3.1. Data Analysis 

Next, we apply our analytic method to three high- 
resolution, high signal-to-noise ratio (S/N) spectra: 
Q 1422+231 at z = 3.62 (Rauch et al. 2001); Q 0130-4021 
at z = 3.023 and Q 0741+4741 at z = 3.21 (Kirkman et al. 
2000, 2005). These quasar were observed with the Keck 
High Resolution Echelle Spectrometer (HIRES), at S/N 
of > 50 - 70. The resolution for Q 1422+231 is - 23.3ft- 1 
kpc, and for the other two is ~ 8 km s . We refer read- 
ers to Rauch et al. (2001) and Kirkman et al. (2005) for 
details of continuum fitting. By fit artificial spectra for 24 
QSOs, Kirkman et al. (2005) measured a typical error in 
continuum level of 1.2%. The actually error in our sample 
should be much smaller due to the higher S/N . 

We select a portion of the spectrum around z ~ 3 to test 
our analytic model. Specifically, we consider the observed 
wavelength band of 4500 < A < 4800 A, corresponding to a 
comoving path length of ~ 180 ft -1 Mpc. For Q 0130-4021 
and Q 0741+4741, only Lya absorption presents in this 
wavelength range, while for Q 1422+231, Ly/3 will cover 
part of this range (but not Ly7) and will have quite a bit 
more aborption. We split the spectrum into segments with 
a size of 60 h~ Y Mpc each so that we have a total of 9 seg- 
ments for the entire sample. We then normalize each of the 
9 segments separately. We did this by varying the optical 
depth t so that the mean flux F — 0.684, consistent with 
the observed value at z = 3 (McDonald et al. 2000; Kirk- 
man et al. 2005). We did not vary the normalization with 
wavelength. We measure gaps in each segment, and the 
errors are estimated by using the variance between these 
9 subsamples. 

We directly measure dP/dl{> L) from the observed 
spectrum. The flat absorption-free gaps are selected by 
searching for segments with transmission F > F t h where 
F t h = 0.99 is our fiducial threshold. Other choices of the 
threshold transmission may be considered plausible, since 
Adelberger et al. (2003) found that the mean transmitted 
flux is > 0.8 within 0.5 ft -1 comoving Mpc of the LBGs and 
possibly higher close to the LBGs. Note that since LBGs 
are selected by having ongoing starburst activity, their ac- 
tive wind may add some absorption by outflowing material 
inside their bubbles. Other galaxies are most likely to be 
surrounded by empty relic bubbles bounded by a thin shell 
of swept-up gas, long after their powerful starburst activ- 
ity has ended. However, a low value of F t h is very likely to 
mix the outflow-driven bubbles discussed here with high- 
transmission gaps caused by gravitationally-induced voids 
(which are undcrdense but not completely empty). We 
have chosen the threshold of 0.99 based on a detailed 
analysis of the results of a numerical simulation from a 
Particle-Mesh (PM) code involves gravity only (see. e.g., 
Meiksin & White 2003; White 1999). We found that the 
PM simulation provides the same distribution of gaps as 
the data up to a threshold of ^0.97, and the data shows 
a clear over-abundance of gaps at a threshold of 0.99: at 
F t h — 0.99 the gravity-induced gaps accounts for less than 
~ 10% of the total gaps found in the observation data. 

In principle we could have used the raw spectrum, the 
noise in each pixel generates artificial gaps with a width 
of just a few pixels. Instead, we use the fitted spectrum, 
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which smoothes out these random pixel noise. Using high- 
order statistic method, Fang & White (2004) found that 
the coherence scale length of spatial correlation in Keck 
data must be smaller than 0.3/i _1 Mpc, which roughly 
corresponds to 12 pixels in the spectrum. So the raw spec- 
trum is smoothed by the boxcar method with a width of 
12 pixels to ensure that we can smooth Poisson noise in 
each pixel but preserve the large scale features. The 1-er 
error bars are obtained directly from Poisson statistics. 

We fit the observed distribution of gap sizes using equa- 
tion (3), with two free parameters: M m i n and Rb- We use 
a least-x 2 analysis with 

X 2 =J2[y(M min ,R b )- yi } 2 /al (4) 

i 

Here y(M m i n , R) is the expected value of dP/dl based on 
equation (3), yi is the observed value of dP/dl, and o, 
is the 1-cr error for each data point. The x 2 fit yields 
the best-fit parameters of M m j n = 1.2 x 10 10 M Q and 
Rb = 480 h~ x comoving kpc, with a x 2 value of 1.2 per 
degree-of-freedom. The model fits the observation data 
well between 0.1 and 1 h~ x Mpc 4 , as indicated by the red 
solid line in Figure 1. 
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Fig. 1. — Cumulative distribution of absorption-free gaps with 
comoving length > L per comoving h^ 1 Mpc. The three heavy 
lines show results from our a 0.48 h^ 1 Mpc and different value of 
M m i„. Data is shown as points (extracted with F t ^ = 0.99) with 
I-ct Poisson error bars. The middle dashed is our best-fit model. 

3.2. Results 

The best-fit value for the minimum halo mass that is ca- 
pable of hosting outflows (due to star formation or quasar 
activity), is consistent with the value expected from hy- 
drodynamic simulations and semi- analytic analysis. The 
mass of a dark matter halo is related to its circular velocity 
V C i rc and collapse redshift z c (3> 1), 

The infall of gas onto dwarf galaxies is significantly sup- 
pressed through photoionization heating of the IGM by 
the cosmic ultraviolet background radiation after cosmo- 
logical rcionization (Rees 1992; Efstathiou 1992; Thoul & 
Weinberg 1996). As soon as the IGM is heated to a tem- 
perature of 1-2 x 10 4 K, its thermal pressure suppresses 
gas infall into shallow gravitational potential wells. The 
minimum halo mass due to photo-ionization heating cor- 
responds to a virial temperature of ~ 10 5 K. The reason is 

4 The correlation function of galaxy halos of ~ 10 10 Mq at z ~ 3 impli 
kpc around these galaxies, is small in our data set. 



that photo-ionization heats the IGM to ~ 10 4 K but when 
the gas virializes into a halo, its density increases by a fac- 
tor of ~ 200 and so its temperature rises adiabatically by 
a factor ~ 200 2 / 3 ~ 35. Only if this final temperature is 
below the virial temperature, will the gas condense by this 
factor of ~ 200. The minimum halo mass for galaxies is 
defined as the threshold where ~ 50% of the baryons that 
were supposed to condense made it. This gives a mini- 
mum virial temperature of ~ 10 5 K. The corresponding 
minimum circular velocity of the halo is (Barkana & Loeb 
2001), 

V mm = 60(T/10 5 K)Hms- 1 . (6) 

Detailed hydrodynamic simulations (e.g., Thoul & Wein- 
berg 1995, 1996; Navarro & Steinmetz 1997; Dijkstra et al. 
2004) confirm this value, showing that at redshift z = 3 
only ~ 50% of the available gas accretes onto halos with 
Vdrc ~ 50kms _1 , and the infall of gas onto halos with 
V C i rc < 25-30 kms" 1 is completely suppressed. This yields 
a minimum halos mass of ~ 10 10 Mq based on equa- 
tion (5), in agreement with the value derived from the 
statistics of absorption-free gaps in the Lya forest. 

The typical bubble radius Rb we derived is consistent 
with both semi-analytic analysis and other observations. 
Based on a starburst wind model, Furlanetto & Loeb 
(2003) showed that at redshift z ~ 3 and halos with masses 
between - 3 x 10 8 M and 10 10 M , the physical (proper) 
radius R b /(1 + z) ranges between ~ 50 to ~ 100 kpc (see 
the upper left panel of their Fig. 1). For halo masses above 
10 10 Mq, the maximum radius can be achieved by the wind 
is nearly constant, ~ 100 kpc. Observations indicate the 
lack of H I within a radius of 500 comoving kpc around 
LBGs, which translates to ~ 125 proper kpc at z = 3 
(Adelberger et al. 2003), although the interpretation may 
be complicated (Kollmeier et al. 2003). 

4. DISCUSSION 

The gas that is swept-up from the bubble volume by 
a galactic outflow, is expected to pile-up in a thin shell 
on the bubble wall, and introduce enhanced absorption at 
both ends of its absorption-free segment in the Lya for- 
est (Furlanetto & Loeb 2003). Indeed, the spectrum of 
Q 1422+231 indicates that almost all of the flat trans- 
mission segments have deep absorption features at one or 
both of their boundaries. Figure 2 shows four examples 
of absorption-free segments with enhanced absorption at 
both ends. In all panels, the upper dotted horizontal line 
show the flux threshold F t h — 0.99, and the lower horizon- 
tal line indicates the mean transmitted flux of F = 0.684 at 
z ~ 3 (bottom). We note that in analogy with supernova 
remnants (Chevalier, Blondin, & Emmering 1992), the 
shell of the swept-up material may fragment into clumps 
due to convective or Raylcigh- Taylor instability. In such 
flat (absorption-free) segment would be bounded 
by a narrow, deep absorption feature only if the line-of- 
sight happens to cross a clump of gas. We find that out of 
23 absorption-free segments with L > 0.1/i _1 Mpc, 6 are 
bounded at both ends by deep absorption features extend- 
ing under F, 12 have such a feature at one end, and only 

that the chance for overlap of bubbles with a radius 480 h~ x comoving 
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5 are not bounded at all by a flux depression under F. 
There is tentative evidence for dumpiness of the bubble 
walls, although the last class of segments may also receive 
a contribution from ionization (Stromgren) spheres around 
bright galaxies. 
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Fig. 2. — Four examples of absorption-free segments in the spec- 
tral interval between 4740Aand 4820Aof Q 1422+231. The raw 
data is indicated by the black fluctuating lines and the smoothed 
spectra are delineated by the red smooth curves. In each panel, 
the top horizontal dotted line marks the threshold F th = 0.99 and 
the bottom horizontal line indicates the mean transmitted flux at 
z ~ 3. All four examples show enhanced absorption at the two edges 
of the central flat (absorption-free) segment, as expected from the 
pile-up of hydrogen on the surface of the associated bubbles. Some 
absorption-free segments (not shown here) are not bounded by deep 
absorption features at both ends, indicating dumpiness of the bub- 
ble material. 

Our analysis assumed that the bubble size can be trans- 
lated to the spectral width of the absorption-free gaps for a 
pure Hubble flow, and ignored the possible effects peculiar 
velocities. In reality, the spectral width may be increased 
due to the bubble expansion relative to the local Hubble 
flow. This would tend to increase the inferred value of Rb 
relative to its true physical size. The local hubble flow at 
z ~ 3 is Vh ~ 32 kms" 1 for bubbles with a radius of ~ 100 



proper kpc. The upper left panel of Figure 2 in Furlan- 
etto & Loeb (2003) shows that for most halos with halo 
mass between 10 9 -10 10 M the average bubble expansion 
velocity is around v exp ~ 40kms _1 . Since the peculiar 
velocity, (v exp — Vh), is much smaller than the hubble flow 
Vh , we can safely ignore its contribution to the broadening 
of the flat spectral segments. On the other hand, neutral 
hydrogen outside the bubble may also have effect in trans- 
lating the bubble size into the spectral width of the gap, 
depending on whether the bubble is expanding into an in- 
fall region (surrounding an overdense environment) or an 
outflow region (in a locally underdense region): the infall 
can require the true physical bubble size to be larger (see, 
e.g., Kollmeier et al. 2003). 

In this Letter we have only attempted to demonstrate 
the feasibility of our method at z ~ 3. The statistical 
significance of our results can be enhanced through the 
analysis of additional data sets of high-resolution quasar 
spectra. We also need better calibrated data, with accu- 
rate continuum level, and more realistic simulatoins. In 
the future, our method can be used to measure the varia- 
tion of M m ; n and Rb with redshift, and to gauge the im- 
pact of galaxy feedback on estimates of the dark matter 
power-spectrum on small scales. 
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